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Abstract 

The equivalence of two complete sets of Poisson commuting Hamiltonians of the 
(super)integrable rational BC n Ruijsenaars-Schneider-van Diejen system is established. 
Specifically, the commuting Hamiltonians constructed by van Diejen are shown to be 
linear combinations of the Hamiltonians generated by the characteristic polynomial of 
the Lax matrix obtained recently by Pusztai, and the explicit formula of this invertible 
linear transformation is found. 
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1 Introduction 


Integrable many-body systems of Calogero-Moser type crop up in a wide range of physical 
applications and are intimately related to important fields of mathematics m ca ei a e]. 
They occur in rational, trigonometric/hyperbolic and elliptic families according to the func¬ 
tional form of the Hamiltonian that inherently also involves a crystallographic root system. 
A further significant feature is the existence of interesting deformations and extensions main¬ 
taining integrability, as is exemplified by relativistic |6j and spin Calogero-Moser systems [7]. 
In this paper we shall deal with the BC„ generalization of the relativistic rational Ruijsenaars- 
Schneider system, which is the simplest member of the systems discovered by van Diejen M. 
We shall stay at the level of classical mechanics, where the BC r) rational ‘RSvD system’ is 
defined by the HaniiltoniarQ 
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( 1 . 1 ) 


Here p, u, k are real parameters for which we impose the conditions p ^ 0, u ^ 0 and 
uk > 0. The generalized momenta 9 = (9±,..., 9 n ) run over M n and the ‘particle positions’ 
A = (Ai,..., A n ) vary in the Weyl chamber 

c = {x G R n | X\ > ■ ■ ■ > x n > 0}. (1.2) 

In the work 0HO] hrst the commutativity of n quantum Hamiltonian difference operators 
was proved. It was then shown mm that the classical limit yields a Poisson commuting 
family having the right functional rank for a Liouville integrable system. Except for the 
rational case, it is still an open problem to generate the classical Hamiltonians of van Diejen 
from a Lax matrix, which would provide a useful tool for analyzing the dynamics of these 
systems. A Lax matrix whose trace is the rational RSvD Hamiltonian (II.ip and whose higher 
spectral invariants provide n independent commuting Hamiltonians was recently found by 
Pusztai H2j. In the papers m S3! the action-angle duality between the hyperbolic BC n 
Sutherland system and the rational BC„ RSvD system was also explored together with the 
scattering properties of these systems. 

The question we answer in this paper is the following. What is the relationship between 
the commuting Hamiltonians introduced by van Diejen and the ones generated by Pusztai’s 
Lax matrix? Both commuting families contain the ‘main Hamiltonian ’ m and exhibit 
rational dependence on the positions and exponential dependence on the momenta. Thus 
one strongly expects that these two sets of commuting Hamiltonians can be expressed in terms 
of each other. Nevertheless, the question appears to be non-trivial since the Hamiltonian H 
is maximally superintegrable na, which entails that it is the member of several inequivalent 
families of n functionally independent functions in involution. 

Here, we shall demonstrate that the n Hamiltonians of van Diejen are linear combinations 
of the coefficients of the characteristic polynomial of the Lax matrix of [12] • The transfor¬ 
mation between the two sets will be shown to be invertible, and its explicit form will be 

deformation parameter /3 > 0 can be introduced by setting H / 3(X 1 9) := (38). Taking Taylor 

expansion of Hp in (3, the leading term reproduces the usual BC„ rational Calogero-Moser Hamiltonian. 
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given as well. Our arguments will rely on the action-angle map constructed with the help of 
Hamiltonian reduction in [12]. This is captured by a symplectomorphism 


S: c x IR n —> c x 


( n \ n 

E dX k A d9 k = E dq k A dp k , (1.3) 

fc=l ' k= 1 


such that H o S depends only the action variables q k . (The notation fits the fact that the 
components of q serve as position variables for the dual system.) Although an explicit formula 
of the action-angle map is not available, we can compute the action-angle transform of the 
commuting Hamiltonians of interest by utilizing that [13] the //-trajectory (A (t),9(t)) with 
initial condition (A, 6) € c x M" has the t —¥ oo asymptotics 


A k {t) ~ t sinh (q k ) - p k and 9 k (t) ~ q k , k = 1,..., n, 


(1.4) 


with ( q,p ) = 5 _1 (A, 6). This will allow us to eventually show that the two sets of Hamilto¬ 
nians at issue correspond to two generating sets of the Weyl group invariant polynomials in 
the variables e ±qk (restricted to the Weyl chamber c). As was already mentioned, we shall 
also find the explicit relationship. As a byproduct, we obtain an algebraic formula for the 
characteristic polynomial of Pusztai’s Lax matrix, which generalizes well-known determinant 
identities for Cauchy-like matrices. 

Section [2] describes the two families of Hamiltonians in play and specifies how they share 
the main Hamiltonian H (11.11) . Our contribution is given by Proposition 1, Proposition 2 and 
Remark 3 in Section [3] Section [4] offers a short discussion of the results and open problems. 
There is also an appendix, where a useful formula of [9] is presented. 

2 Two families of commuting Hamiltonians 


2.1 Hamiltonians due to van Diejen 

In |8[ [TT] the following complete set of Poisson commuting Hamiltonians was given: 
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( 2 . 2 ) 


It is worth noting that J c in (12.ip denotes the complementary set, and the contribution to 
Hi coming from J = 0 is U^ i. The relatively simple form of Uk, p above was found in HU- 
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Equation (12. ip makes sense for l = 0, as well, giving H 0 = 1. In the rational case the functions 
v and w take the following forixH 


v(x) 


x + ifi 
x 


w(x) 


x + w 

X + IK 

X 

X 


(2.3) 


Up to irrelevant constants, Hi reproduces the Hamiltonian H (11. lj) . Indeed, one can check 
that Hi = 2(H — n). 

Take any (A, 6) e c x M n , set (q,p) = S^ 1 (X,6) and consider the //-trajectory (A (t),9(t)) 
with initial condition (A ,9). Notice that the Hamiltonian Hi (12.ip is constant along the 
//-trajectory. By utilizing the asymptotics (ll.4j) . one can readily check that 


ff,)(?,p)= limi/,(A((), «(*))= V (-2)‘-M 

t —^oo z—/ 

JC{l,...,n}, |J|<( 

Sj=±l, j&J 


From now on we let Hi stand for the pullback S*Hi just computed, and stress that it depends 
only on the variable q. 



("-p, l ) C08t(9 "' ) - (24) 


2.2 Hamiltonians obtained from the Lax matrix 

We recall some relevant objects of D2t First, prepare the 2 n x 2 n Hermitian, unitary matrix 


C = 


On In 
In 0 n 


and the 2 n x 2 n Hermitian matrix 

h( A) = 


a(diag(A)) /(diag(A)) 
—/(diag(A)) a(diag(A)) 


containing the smooth functions a(x), b(x) given on the interval (0, oo) C I by 
f , Vx+ Vx 2 + K 2 .... 1 1 

a[x) = --, b[x) = i k- 


y/2x \Jx + \JX 2 + K, 2 
Then introduce the vectors z(A) € C n , F(\,6) e C 2n by the formulae 


z i(ty ~ ~ 
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and 


F,(A,0) =e-2|z t (A)|2, F„ +i (A,0) =^(A)F z (A,0)" 1 , 

l — 1,..., n. With these notations at hand, the 2 n x 2n matrix 

■W) = 

1 h + — Afc 


(2.5) 


( 2 . 6 ) 


(2.7) 


( 2 . 8 ) 


(2.9) 


( 2 . 10 ) 


2 The parameters appearing in 8j[TI] can be recovered by introducing /3 as in footnote 1 and then writing 
/r, Ho, Ho l° r Pn, fiv, fin, respectively. In the convention of H3j, our n, 9 and q correspond to 2/x, 29 and 2 q. 
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( 2 . 11 ) 


with A = diag(A, —A) is used to define the ‘RSvD Lax matrix’ |12|: 

L(\,6) = h(A)- 1 A(A,0)h(A)" 1 . 

The matrices h, A, and L are invertible and satisfy the relations 

ChC = h~\ CAC = A~\ CLC = L~ l . (2.12) 

Their determinants are 

det(h) = det(A) = det(A) = 1. (2.13) 

Let K m denote the coefficients of the characteristic polynomial of L (12.lip . 

det(A(A, 9 ) — xl 2 n ) = Lq(A, 9)x 2n + AA(A, 9)x 2n 1 + • • • + K 2 n — i(A, 9)x + A 2 n (A, 9). (2.14) 

An immediate consequence of (l2.12IUI2.13p is that 

-^ 2 n—m = A );m 0,1,..., n, (2.15) 

thus the functions K 0 = l,Ki,..., K n fully determine the characteristic polynomial (I2.14p . 
The first non-constant member of this family is proportional to H (II.ip . that is K\ = —2 H. 
The asymptotic form of the Lax matrix L (12.131) is the diagonal matrix 

diag(e _<? , e q ), (2.16) 


hence the action-angle transforms of the functions K m (m = 0,1,... ,n) can be easily com¬ 
puted to be 


(S*K m )(q,p) = (-iy 


I — I 
L 2 J 

E 


E 


a=0 Jc{l,...,n}, \J\=m—2a 
£j=±i, jeJ 


n~\J\ 


cosh(g, 


eJ ) 


(2.17) 


Of course, we used the asymptotics (11.41) and that K m is constant along the flow of H. Now 
we introduce the shorthand /C m := S*K m , and observe that it only depends on q. 


3 Relation between the two families of Hamiltonians 

It is worth emphasizing that finding a formula relating the families {H{\\ L 0 and {K m }^ n=0 is 
equivalent to finding a relation between their action-angle transforms {'Hi}} L 0 and {/C m }^ =0 . 

Proposition 1. There exists an invertible linear relation between the two families {Hi}\ L 0 
and 

Proof. Let us introduce the auxiliary functions 

Mk{q) — cosh(g £ j), q G M n , k = 0,1,..., n. (3.1) 

JC{l,...,n}, \J\=k 
£j=±l, j£j 

For any l G {0, 1,... , n} the Hamiltonian Hi (12.41) is a linear combination of Aio, A4i,. . ., Aii, 

k 
k 
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n — 
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This shows that the matrix of the linear map transforming {A4k}k =o i n t° {^- 1 }?= o is lower 
triangular with ones on the diagonal, hence the above relation is invertible. Similarly, 
any function JC m (j2.17j) . m G {0,1,..., n) can be expressed as a linear combination of 
A4 m , Ai m - 2 , • • •, or A4 m , Ai m - 2 ,..., AI 2 , A4 0 depending on the parity of m, that is 


LfJ 

IC m (q) = (-l) m Y, 


a=0 


n — (m — 2a) 


A4 m -2 a (q)- 


(3.3) 


Hence the linear transformation relating {A ik} k =o 1° 0 l ias a l° wer triangular matrix 

with diagonal components ±1, implying that it is invertible. This proves the existence of an 
invertible linear relation between the two families {Hi }^ 0 and {/C m }() l=0 . □ 

Now, we prove an explicit formula expressing Hi as linear combination of {JC m } l m=0 . 

Proposition 2. For any fixed n G N, l G {1, ..., n} and q G M n we have 

i -i 


[-!)'%(<?) = r,( q ) + 


2 (n — m ) 


2{n — m) — (/ — m) 


m =0 


(n — l) + in — m) 
l — m 




(3.4) 


Proof. Substitute IC m (I2.17P into the right-hand side of the expression above to obtain 


LfJ 

EE E M) 

k =0 a=0 \J\=k—2a 

£j=±l, jeJ 


k 2 [n - k) f(n — l) + (n - k) 


x 


n — (k — 2a) 


cos. 


2(n — k) — (l — k) 

[L\ 

HQsj) + Y Y 


l - k 


x 


(- 1 ) 


! I n — (l — 2a) 


a=0 \J\=l—2a 

£j= ±1, jeJ 


cosh (q eJ ). (3.5) 


Since k = \ J\ + 2a it is obvious that (—l) fe = (—1) Multiply (13. 5 p by (—l) 1 and change 
the order of summations over a and J to get 


LhrJ 


2[n — (| J| + 2a)] 


V (-l)'-M V 

h 2 l n-(\J\ + 2a))- V -(\J\+2a)] 

ej=± i, jeJ 


X 


X 


(n — l) + (n — (| J\ + 2a))^ fn - \ J\ 
l — (\J\ + 2a) 


cos. 


HQsj) + Y c °sh (q eJ ). (3.6) 


J C{l,...,n}, \J\=l 
£j=±i, jeJ 


Now, comparison of (13. 2 p with (13. 6 p leads to a relation equivalent to (13.41) . 


2[n-(|J|+2o)] 


L^J 

H, (I A + 2a)] - [i - (M + 2a)] 

'2n — {l + | J\ + 2 a)\ fn — \J\\ / fn — |J| 


■x 


l-{\J\ + 2a) 


l-\J\ 


= 2 HJ| . (3.7) 
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For n — l in (13.71) one obtains 



if / — | J\ is odd, 


2 hj| , if l - | J\ is even, 


(3.8) 


which are well-known identities for the binomial coefficients. This means that ( 13 . 41 ) holds for 
l = n for all n G N, which implies that if we consider n + 1 variables it is sufficient to check 
the cases l < n + 1. With that in mind let us progress by induction on n and suppose that 
( 13 . 41 ) is verified for all 1 < l < n for some n G N. 

First, notice that the Hamiltonians Hi (12.41) satisfy the following recursion 


q n , q n+ i) = %i{qu ■ ■ •, q n ) + 4 sinh 2 (^yi)?Zz_i(gi,..., q n ). (3.9) 

This can be checked either directly or by utilizing that Hi is the Z-th elementary symmetric 
function with variables sinh 2 ^) (see Appendix A). Similarly, the functions /C*. (I2.17|) satisfy 


ICkiqi, ■ ■ -AnHn+i) = fck(qi, ■ • .,g n )-2cosh(g n+ i)/C fc _i(gi,..., g„)+/C fc _ 2 (gi, • ■ .,q n ), (3.10) 


with /C_i = 0. Let us introduce some shorthand notation, such as the M /+1 vectors 

H(n) :^(H 0 ,-H 1 ,...,(-l) l H l ) T and K{n) := (ZC 0 , JC U ..., /C,) T (3.11) 

and the R0 +1 ) x ( ,+1 ) matrices 

( 2 (n - k) f{n- j) + (n - fc)\ 

^(n) i+lifc+1 :=h(n-k)-(j-k)\ j-k )' (3.12) 

[o, if j<k, 

where j, k G {0,..., /} and 

H(n, 7i + 1) := 1; +1 -4sinh 2 (^yi)X_ 1 , /C(n,n + 1) := l i+ i - 2 cosh(g n+ i)X_i+X_ 2 (3.13) 

with (X_ m )j + i : fc + i := m > 0. The relations (13.9)) and (13.101) can be written in the 

concise form 


H(n + 1) = H(n, n + l)H(n), JC(n + 1) = JC(n, n + 1)/C(n) (3-14) 

and our assumption is condensed into 

H{n) — A(n)K.(n). (3.15) 

Using this notation it is clear that the desired induction step is equivalent to the matrix 
equation 

H{n , n + l)Xl(n) = A(n + 1 )/C(n, n + 1). (3.16) 
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Spelling this out at some arbitrary (j, fc)-th entry gives us 


A + B A4\ 

a U; 


— 4 sink 2 


(a\ A + B (A+l\ 
\2j A + l V^-lJ 


A + B + 2 (A + 2\ 

A + 2 V B J 


2 cosh (a) 


A + B (A + 1\ A + B -2 


A + 1 \B 


+ 


A 


A 

B-2 


(3.17) 


where 

A := 2n — j — k, B := j — k, a := q n+ \. (3.18) 

A simple direct calculation shows that (13.171) indeed holds implying that (13.41) is also true for 
n + 1 for any l < n. The case l = n + 1 is given by the argument preceding induction. This 
completes the proof. □ 


Remark 3. We showed in Proposition 1 that the relation (13.41) is invertible. Without 
spending space on the proof, we note that the inverse relation can be written explicitly as 



(3.19) 


4 Discussion 

In this paper we demonstrated that the commuting Hamiltonians of the rational RSvD sys¬ 
tem constructed originally by van Diejen are linear combinations of the coefficients of the 
characteristic polynomial of the Lax matrix found recently by Pusztai, and vice versa. The 
derivation utilized the action-angle map and the scattering theory results of dam. Our 
Proposition 2 gives rise to a determinant representation of the somewhat complicated ex¬ 
pressions Hi in (12.11) . It could be of some interest to provide a purely algebraic proof of the 
resulting formula of the characteristic polynomial of the Lax matrix. 

The configuration space c (II.2|) is an open Weyl chamber associated with the Weyl group 
W (BC n ), and after extending this domain all Hamiltonians that we dealt with enjoy W (BC n ) 
invariance. In particular, the sets {'H;}” =0 , {/Q}” =0 and {A4;}jL 0 represent different free 
generating sets of the invariant polynomials in the functions e ±qk ( k = 1 of the 

action variables q\. acted upon by the sign changes and permutations that form ^(BC^). In 
order to verify this, it is useful to point out that the B^BC^) invariant polynomials in the 
variables e ±qk are the same as the ordinary symmetric polynomials in the variables cosh(gfc). 
The statement that {"H/j-jL 0 is a free generating set for these polynomials then follows, for 
example, from the identity presented in Appendix A. 

Analogous statements hold obviously also for the different real form of the complex ra¬ 
tional RSvD system studied in |15], which is also superintegrable. 

An interesting open problem for future work is to extend the considerations reported here 
to the hyperbolic RSvD system having five independent coupling parameters. 
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A Tii as elementary symmetric function 

Fix an arbitrary n G N and l G {0,1,..., n} and let ei stand for the Z-th elementary symmetric 
polynomial in n variables xi ,..., x n , i.e., eo(£i,..., x n ) = 1 and for l > 1 

e-i(x i,...,x n )= ^ x h'" x ii • A- 1 ) 

1 <ji<-<ji<n 

In the text, we referred to the following useful result due to van Diejen ([9] Proposition 2.3). 
For convenience, we present it together with a direct proof. 

Proposition A. By using (12.4p it can be shown that 

Hi{q ) = 4*e z (sinh 2 y,..., sinh 2 y). (A.2) 

Proof. First, ei has the equivalent form 

e/(sinh 2 ..., sinh 2 y) = ^ JJ sinh 2 1-. (A.3) 

JC{l,...,n}, \J\=l j&J 

Utilizing the identity sinh 2 (a/2) = [cosh(a) — l]/2 casts the right-hand side into 

J2 -!] = 2 ~ l II cosh (^)- ( A - 4 ) 

Jc{l,...,n}, |J|=Z jeJ Jc{l,...,n}, |J|=Z Kcj k£K 

The two sums on the right-hand side can be merged into one, but the multiplicity of subsets 
must remain the same. This results in the appearance of a binomial coefficient 


E 


J C{l,...,7l}, |J|<z 


-iy-w\ 

~^2 l 


n-\J\ 

l-\J\ 


JJcosh(^) = 

- E 

JC{l,...,n}, \J\<1 
£j=± 1, j£J 


-1) 


l-\J\ 


\ l 


n 


\J\ 

\J\ 


n cosh(£j5j), (A.5) 
i&J 


where we also used that cosh is an even function and compensated the ‘over-counting’ of 
terms. Now, let us simply pull a A~ l factor out of the sum to get 

A- 1 J2 (— :2 ) HJ l(?_ ij|) II cos h( £ ^i)- ( A - 6 ) 

Jc{i,...,n}, \J\<1 ' ' '' jeJ 

ej=± i, jeJ 


Recall the following identity for the hyperbolic cosine of the sum of a finite number, say A, 
real arguments (see [16] Art. 132 and apply cos(ia) = cosh(a)) 


cosh 



r N 

]J cosh(afc) 

- k=\ 


L?J , 

^ e 2m (tanh(ai),..., tanh(ajv)) , 

771—0 


(A.7) 


where e 2 m are now elementary symmetric functions with arguments tanh(ai),..., tanh(ajv)- 
Note that for any m > 0 and set of signs £ there is another one E, such that 4m = 4m i 
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therefore by using (1A.7|) we see that (1A.6|) equals to 


i-i 


E 




JC{l,...,n}, |J|<Z 

£j=± 1, j&J 


n—\J\ 

l~\J\ 


II cosh(e J g J ) 


j&J 
= 4~ l 


m =0 


E 


s J,£ - 

°2m ~ 




Jc{l,...,n}, |J|<Z 
£j=±i> ieJ 


n- |J| 

l~\J\ 


cosh (g eJ ). (A.8) 


Applying (I2.4|) concludes the proof. 


□ 
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